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Heat  and  mass  transfer  with  electrochemical  reaction  in  an  anode-supported  flat-tube  solid  oxide  fuel 
cell  (FT-SOFC)  is  studied  by  means  of  three-dimensional  numerical  simulation.  The  distributions  of  the 
reaction  fields  in  the  anode-supported  FT-SOFC  are  found  to  be  similar  to  those  in  the  planar  SOFC  with 
co-flow  arrangement.  However,  in  comparison  with  the  latter,  the  concentration  and  activation  overpo¬ 
tentials  of  the  former  can  be  reduced  by  additional  reactant  diffusion  through  the  porous  rib  of  the  fuel 
channel.  Parametric  survey  reveals  that,  for  a  fixed  activation  overpotential  model,  the  output  voltage 
can  be  improved  by  increasing  the  pore  size  of  anode,  while  the  cross-sectional  geometry  has  smaller 
effect  on  the  cell  performance.  Based  on  the  results  of  three-dimensional  simulation,  we  also  develop 
a  simplified  numerical  model  of  anode-supported  FT-SOFC,  which  takes  into  account  the  concentration 
gradients  in  the  thick  anode  of  complex  cross-sectional  geometry.  The  simplified  model  can  sufficiently 
predict  the  output  voltage  as  well  as  the  distributions  of  temperature  and  current  density  with  very  low 
computational  cost.  Thus,  it  can  be  used  as  a  powerful  tool  for  surveying  wide  range  of  anode-supported 
FT-SOFC  design  parameters. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  is  expected  to  be  a  promising  alter¬ 
native  power  generation  system  with  its  high  energy  conversion 
efficiency,  wide  fuel  adaptability,  and  possible  exhaust-heat  uti¬ 
lization  [1].  Recently,  a  flat-tube  SOFC  (FT-SOFC)  has  attracted 
increasing  attention  among  various  cell  configurations  due  to  its 
high  power  density  and  the  relative  easiness  of  fabrication  and  gas 
sealing  [2,3].  For  successful  implementation  of  SOFC,  however,  the 
cost  performance  should  be  further  improved  both  in  the  manu¬ 
facturing  and  operating  processes.  Therefore,  development  of  cells 
with  high  efficiency  and  high  power  density  is  a  current  priority 
issue. 

The  performance  characteristics  of  SOFC  can  be  increased  by 
reduction  of  overpotentials,  which  result  from  electrochemical 
processes  in  the  cell.  One  way  to  reduce  the  overpotential  is  to 
properly  manage  the  heat  and  mass  transfer  associated  with  the 
electrochemical  reaction.  In  order  to  study  such  transport  and  reac¬ 
tion  phenomena,  numerical  simulation  should  be  helpful  to  obtain 
detailed  (and  experimentally  unavailable)  information  in  the  cell. 
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In  the  last  decade,  a  number  of  numerical  studies  have  been 
reported  for  the  detailed  heat  and  mass  transfer  in  SOFCs  [4-7]. 
These  studies  reveal  that  the  flow  and  thermal  fields  in  the  cell  have 
three-dimensional  distribution  strongly  depending  on  the  design 
parameters  of  the  cell  components.  However,  most  of  the  previous 
studies  assume  homogeneous  reaction  field  and  have  not  treated 
the  coupled  transport  and  reaction  phenomena  that  largely  affects 
the  cell  performance  characteristics.  Exceptionally,  several  recent 
investigations  of  tubular,  planar,  and  cathode-supported  FT-SOFCs 
[3,8-10]  account  for  such  interactive  effect. 

When  compared  to  the  cathode-supported  FT-SOFC  mentioned 
above,  the  anode-supported  FT-SOFC  has  potential  advantages  due 
to  its  lower  internal  resistance  and  higher  limit  current  [1,11],  and 
therefore  it  attracts  considerable  attention.  The  anode-supported 
FT-SOFC  has  a  thick  anode,  of  which  cross-sectional  geometry  is 
often  complicated.  Its  performance  is  determined  by  the  three- 
dimensional  heat  and  mass  transport,  which  may  be  strongly 
dependent  on  the  cell  design,  such  as  the  geometrical  configuration 
and  the  microstructure  of  cell  components.  To  the  authors’  knowl¬ 
edge,  however,  detailed  numerical  simulation  of  coupled  transport 
and  reaction  phenomena  in  the  anode-supported  FT-SOFC  has  not 
been  reported. 

The  purpose  of  the  present  study  is  to  investigate  the  funda¬ 
mental  transport  and  reaction  phenomena  in  an  anode-supported 
FT-SOFC.  First,  we  perform  numerical  simulation  with  a  full  3D 
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Nomenclature 

cpj  molar  specific  heat  of  reactant  i  (J  mol-1 I<-1 ) 

Cp  isobaric  specific  heat  (J  kg-1  K-1 ) 
dp  pore  diameter  of  the  electrode  (m) 

Dt  Diffusivity  of  the  reactant  i  (m2  s-1 ) 

E  electromotive  force  (EMF)  (V) 

F  Faraday  constant:  9.649  x  104  (C  mol-1 ) 

AGh2o  Gibbs’  free-energy  change  accompanied  by  water 
formation  (J  mol-1 ) 

h  heat  transfer  coefficient  (W  m-2  K-1 ) 

AHh2o  enthalpy  change  accompanied  by  water  formation 
(J  mol-1 ) 

i  current  density  (A  m-2 ) 

fhj  molar  flux  of  reactant  i  (mol  s)-1 
Mi  molar  weight  of  the  reactant  i  (kg  mol-1 ) 
n  number  of  electrons  participating  in  the  reaction 

Pi  partial  pressure  of  the  reactant  i  (Pa) 

P  pressure  (Pa) 

Q  heat  generation  in  the  unit  volume  of  the  electrodes 

(Wnrr3) 

R0  gas  constant:  8.314  (J  mol-1  K-1 ) 

Ri  electrical  resistance  of  cell  component  i  (£2) 

T  temperature  (I<) 

u  flow  velocity  (ms-1) 

Vceli  output  voltage  of  the  cell  (V) 

Wcell  volumetric  power  density  of  the  cell  ( W  m-3 ) 
x,  y,  z  Cartesian  coordinates 

X/  molar  fraction  of  the  reactant  i 

Yi  mass  fraction  of  the  reactant  i 

Greek  symbols 

e  porosity  of  the  electrode 

?7act  activation  overpotential  (V) 

Peon  concentration  overpotential  (V) 

dohm  Ohmic  loss  (V) 

k  flow  permeability  in  the  electrode  (m2) 

k  thermal  conductivity  (W  m-1  K-1 ) 

fi  viscosity  (Pas) 

p  density  (kg  m-3) 

a  resisivity  (Q  cm) 

r  tortuosity  factor  of  the  electrode 

Subscripts 
an  anode 

avr  average 

bulk  bulk  average  in  the  cross-sectional  area  of  the  cell 

ca  cathode 

ch  fuel  channel 

ele  electrolyte 

gas  gas 

gen  generation  accompanied  by  electrochemical  reac¬ 

tion 

H2  hydrogen 

H20  water  vapor 

IB  interface  boundary  between  solid  and  porous  com¬ 

ponent  of  the  cell 
in  inlet  of  the  cell 

int  interconnector 

02  oxygen 

sol  solid 

TPB  three-phase  boundary 


model  and  discuss  the  detailed  characteristics  of  flow  and  reaction 
fields  in  the  cell.  The  effects  of  the  design  parameters  of  the  cell 
are  surveyed,  especially  focusing  on  the  cross-sectional  configura¬ 
tion  and  the  microstructure  of  the  porous  anode.  We  also  propose 
and  assess  a  simplified  numerical  model  of  the  anode-supported 
FT-SOFC  for  practical  design  of  fuel  cell  systems  [1]  based  on  the 
knowledge  obtained  from  the  3D  simulation. 

2.  Numerical  procedure  for  three-dimensional  simulation 

2.1.  Computational  model 

Fig.  1  shows  the  anode-supported  FT-SOFC  considered  in  the 
present  study.  As  shown  in  Fig.  1(a),  it  has  a  periodic  structure  in  the 
spanwise  (z-)  direction.  The  unit  cell  is  composed  of  an  electrolyte, 
a  cathode,  an  anode  with  a  fuel  channel,  an  interconnector,  and  an 
air  channel  between  the  adjacent  layers  of  cells.  The  thick  anode 
plays  a  role  to  mechanically  support  the  cell.  The  air  and  fuel  flow 
directions  are  in  co-flow  arrangement. 

The  geometry  under  computational  study  is  depicted  in  Fig.  1  (b). 
We  take  only  a  half  of  the  unit  cell  shown  in  Fig.  1(b),  because 
the  flow  in  each  channel  is  laminar  and  plane-symmetric  since 
the  Reynolds  number  is  low  (Re  =  1-50)  under  all  operating  con¬ 
ditions.  The  interface  between  the  electrolyte  and  electrodes  are 


Fig.  1.  (a)  Schematic  of  the  anode-supported  FT-SOFC.  (b)  Computational  model 
of  the  anode-supported  FT-SOFC.  The  thin  layers  of  electrolyte,  cathode,  and  inter¬ 
connector  (of  which  thicknesses  are  20  |xm,  100  p,m  and  100  |xm,  respectively,  and 
treated  as  boundaries  in  the  simulation)  are  drawn  not  in  scale  so  as  to  be  visible. 
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Table  1 

Physical  properties  of  cell  components.  The  values  are  referring  to  the  previous  studies  [1,9,12,13] 
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treated  as  the  three-phase  boundary  (TPB).  The  electrodes  are 
porous  media.  The  electrolyte  and  the  interconnector  are  treated  as 
solid  boundaries,  which  couple  the  fuel  and  air-sides  as  described 
below.  The  parameters  of  geometrical  configuration  to  be  changed 
in  the  present  study  are  the  rib  width,  Zzrib,  the  rib  thickness,  ly>an, 
whereas  the  width  and  height  of  the  fuel  channel,  Zz  ch  and  ly  ch,  are 
fixed. 

The  physical  properties  of  each  component  are  listed  in  Table  1. 
The  resistivity  is  considered  as  a  function  of  temperature,  whereas 
other  properties  are  assumed  to  be  constant  [1,9,12,13].  Among  the 
material  properties  of  anode,  the  porosity,  sa n,  and  the  pore  diame¬ 
ter,  dPian,  are  varied  in  the  present  study.  The  tortuosity  factor  of  the 
porous  media,  r,  is  correlated  with  the  porosity  [14].  Note  that  the 
tortuosity  factor  of  anode,  ran,  used  in  the  previous  studies  of  SOFC 
has  large  scatter  [4,9,11,15-18].  In  the  present  study,  we  assume  a 
linear  function  as  follows: 

Tan  =  5.0  —  4£an-  (1) 

For  the  cathode,  the  porosity  and  the  tortuosity  factor  (£ca  and 
rca)  are  fixed  as  £ca  =  0.5  and  rca  =  3.0,  respectively. 

The  electrochemical  reaction  is  calculated  by  using  an  equiv¬ 
alent  circuit  as  shown  in  Fig.  2.  This  model  consists  of  the 
electromotive  force  (EMF)  and  activation  overpotential  generated 
in  the  cell,  internal  resistance  of  the  cell  components  and  the 
current  path.  The  cathode  and  the  interconnector  are  treated  as 
equi-potential  surfaces.  The  average  current  density  of  the  elec¬ 
trolyte,  ieie.ave*  which  is  defined  as  the  current  density  averaged 
over  the  whole  area  of  electrolyte,  is  kept  constant,  while  the  local 
current  density  varies  in  space  and  is  determined  as  a  result  of  the 
heat  and  mass  transfer  as  well  as  the  electrochemical  reaction.  Since 
the  steady-state  performance  of  the  cell  is  discussed  in  the  present 
study,  the  capacitance  between  the  electrolyte  and  the  electrodes 
is  neglected. 


2.2.  Governing  equations  and  boundary  conditions 

2.2.1.  Momentum  transport 

The  transport  equations  for  mass  and  momentum  are 

+  V  •  (epu)  =  0,  (2) 

^£^1  +  v  ■  {(spu)u)  =  -v  (sP  +  •  u) 

2 

+  V  •  [e/L  {vu  +  T(Vu)}]  -  ^-^u.  (3) 

These  are  the  volume-averaged  conservation  equations  of  mass  and 
momentum  in  the  porous  media  (see,  e.g.,  [19]).  Flere,  Darcy’s  law 
is  adopted  to  describe  the  superficial  flow  velocity  in  the  porous 
media  (i.e.,  the  last  term  in  Eq.  (3)). 

Fully  developed  laminar  velocity  profiles  are  given  at  the  inlets 
of  air  and  fuel  channels.  The  flow  rates  are  determined  by  the  pre¬ 
scribed  air  and  fuel  utilization  rates  shown  in  Table  2.  The  inlet 
gas  pressure,  Pgas,in*  is  fixed  at  1  atm.  Zero-flux  boundary  condition 
is  applied  at  the  symmetric  plane.  No-slip  boundary  conditions  are 
applied  on  the  impermeable  boundaries  such  as  the  electrolyte  and 
the  interconnector. 


2.2.2.  Transport  of  chemical  species 

The  transport  of  chemical  species  is  governed  by  the  following 
equation: 


+  V  ■  ( spuYi )  =  V  .  (spDj  gffVYj).  (4) 

Here,  the  effective  diffusivity  of  reactant  i  in  the  porous  media 
(Dieff)  is  calculated  by  using  the  dusty-gas  model  [4,20],  which 
reads 


^i,  eff  — 


e  A  -  V,(1  -  —) 


-1 


Df 


Di 


i,  k 


(5) 


where  Di  m  and  Di  k  are  the  multi-component  molecular  diffusivity 
and  the  Knudsen  diffusivity  of  reactant  i,  respectively,  while  Mmix 
represents  the  averaged  molar  weight  of  the  reactant.  Note  that  Dim 
and  Dj  k  are  functions  of  the  microstructure  of  the  electrode,  local 
temperature,  T,  and  the  reactant  concentration,  Yz-  as  well  as  the 
other  physical  properties  of  the  fluid  [4,15,16,21].  The  summation 
of  Eq.  (4)  for  all  reactants  is  identical  to  Eq.  (2). 

Methane  completely  reformed  at  a  steam  to  carbon  ratio  of 
2.5  is  selected  as  the  inlet  composition  of  fuel.  Namely,  the 
molar  fractions  are  XH2)in  =  8/11,  ^H2o,in  =  1/11,  and  ^co2,in  = 


Table  2 

Parameters  common  to  all  the  cases 


Inlet  gas  temperature,  Tgas  in  (°C) 

600 

Inlets  pressure,  Pgas,in  (atm) 

1.0 

02  Utilization  rate 

0.3 

H2  Utilization  rate 

0.85 

Width/height  of  fuel  channel  (mm) 

1.0 

Cell  length  (mm) 

100 
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2/11,  respectively.  The  corresponding  mass  fraction  is  given  at  the 
inlet.  The  boundary  condition  on  the  TPB  is  given  by 


with  the  values  of  resistivities,  cr,  shown  in  Table  1.  The  electromo¬ 
tive  force  (EMF),  E,  is  computed  by  the  Nernst  equation,  i.e., 


d(spYi)  _  _j_  ie\eMf 

9y  TPB  ~~  nFDUeff 


IE  O2,  H2,  H2O, 


(6) 


where  n  =  4  for  O2  and  n  =  2  for  H2  and  H20.  The  sign  of  the  right- 
hand-side  term  is  plus  for  02  and  H20,  and  minus  for  H2.  The 
concentration  gradient  of  N2  and  C02  on  the  TPB  is  assumed  to 
be  zero,  since  these  gases  are  inert.  On  the  interconnector  sur¬ 
face,  which  is  impermeable,  the  concentration  gradient  of  all  the 
reactant  should  be  zero. 


2.2.3.  Heat  transport 

The  energy  transport  equation  for  a  porous  media  as  well  as  for 
the  air  and  fuel  channels  is  given  by  using  the  local  equilibrium 
temperature  as: 


d[{£PgasCp  ,gas  +  (l  £)Psol^p,sol}^]  ,  _  ,  r 
- — - b  V  •  (,£pgasCp,gas 

=  V  •  [{£Ap,gas  +  (1  -  £)Ap?sol}VT]  +  Q. 


UT) 


(7) 


This  is  a  volume-averaged  equation  of  the  heat  conduction  in  the 
solid  part  and  the  energy  transport  in  the  fluid  phase.  The  temper¬ 
ature  difference  between  the  two  phases  in  the  porous  electrode 
can  be  neglected  under  the  present  condition.  The  local  heat  pro¬ 
duction,  Q,  in  Eq.  (7)  is  determined  by  the  local  ohmic  loss  in  the 
electrode. 

The  solid  and  porous  boundaries  at  the  cell  inlet  and  outlet 
are  assumed  to  be  adiabatic.  The  inlet  gas  temperature,  Tgas  in,  is 
600  °C  throughout  the  present  study  unless  specified  otherwise. 
We  neglect  the  temperature  distribution  inside  the  thin  compo¬ 
nents  (i.e.,  the  electrolyte  and  the  interconnector)  in  the  direction 
of  thickness,  because  they  are  much  thinner  than  the  electrodes 
and  the  fluid  passages.  Namely,  the  temperatures  on  both  sides, 
TIB  air  and  TIB  fuel  (where  the  subscript  IB  denotes  the  boundary  of 
the  electrolyte  or  interconnector  and  air  and  fuel  denote  the  air 
and  fuel  sides,  respectively)  and  that  in  the  thin  component,  Tsoh 
are  assumed  to  be  the  same,  i.e., 


Tib, air  =  Tib, fuel  =  Esob  (8) 

and  it  is  computed  by  solving  the  following  set  of  equations  together 
with  Eq.  (8): 


E 


2F  4F  \  Ph20  j 


(11) 


The  partial  pressure  of  species  2,  p,-,  is  obtained  as  a  product  of  the 
total  local  pressure  and  the  molar  fraction,  Xif  where  the  molar 
fraction  is  calculated  by  using  the  mass  fraction,  Yit  and  the  molar 
weight,  Mj,  as 


Yj/Mi 


(12) 


In  the  present  study,  the  concentration  overpotential  is  defined 
as  the  difference  between  the  EMF  computed  based  on  the  mean 
concentration,  Ebuik,  and  the  EMF  based  on  the  local  concentration 
on  the  TPB,  Etpb>  be., 


Peon  =  Ebuik  -  Etpb 


f  (Po2 , bulk/ Po2  ,TPB  )1  /2 (Ph2  ,bulk/PH2  ,TPB  ) 
2E  j  (PH20,bulk/  Ph20,TPb) 


(13) 


The  activation  overpotentials  in  the  cathode  and  anode  are 
implicitly  given  by  the  Butler-Volmer  equation,  i.e., 

where  a  is  the  transfer  coefficient  (0  <  a  <  1 )  and  n  is  the  number  of 
electrons  participating  in  the  reaction  (n  =  4  for  cathode  and  n  =  2 
for  anode),  respectively.  In  the  present  study,  a  is  assumed  to  be  0.5. 
The  exchange  current  densities,  i0,  are  expressed  as  follows  [22] 


where  y  and  Eact  are  the  pre-exponential  factor  and  the  activation 
energy  of  the  electrode,  respectively.  In  the  present  study,  the  values 
of  these  parameters  are  chosen  as  shown  in  Table  3  following  the 
previous  studies  [9,13,22]. 

Finally,  the  output  voltage  of  the  cell,  Vcen,  can  be  obtained  by 


,,  dT 
Pair  ~  ±yMB,air 


Pfuel  —  T^IB,fuel  ^ 
,  Qair  +  Pfuel  =  Pgene 


IB,  air 

dT 


IB,  fuel 


0) 


Here,  AIB  is  the  thermal  conductivity  of  fluid,  while  pair  and  qfuei 
denote  the  fractions  of  qgen  distributed  to  the  air  and  fuel  sides, 
respectively.  The  heat  generation  in  the  unit  area  of  the  electrolyte 
and  the  interconnector,  qgen,  is  given  by 


Q  gen.ele— *ele  \  2  p 
Pgen,int  =  *elePohm,int’ 


AH, 


h2o 


-E 


+  Pohm,ele  +  Pact.ca  +  Pact, an] 


,  (10 


and  the  division  of  qgen  into  qair  and  qfuel  follows  the  local  temper¬ 
ature  gradient  determined  through  the  iteration  between  the  heat 
and  mass  transports  and  the  electrochemical  reaction. 


2.2.4.  Electrochemical  reaction 

The  local  current  density,  i,  and  the  ohmic  losses,  pohm,  of  the 
cell  components  are  computed  by  using  Ohm’s  and  Kirchhoff  s  laws 


VCell  —  Ebulk  —  (Peon  +  Pohm  +  Pact)-  (16) 

2.3.  Numerical  procedure 

Eqs.  (2)-(4)  and  (6)  are  solved  by  developing  a  simulation 
code  based  on  the  standard  SMAC  procedure  [23]  under  the 
low-Mach  number  approximation  [24].  Namely,  the  second-order 
accurate  central  difference  method  is  used  for  the  spatial  dis¬ 
cretization  on  the  staggered  mesh  system  and  time  integration 
is  done  by  using  the  delta-form  fractional  time  step  method. 
The  computational  grid  has  a  constant  spacing  in  each  direction: 
Ax  =1.5625  mm,  Ay  =  0.1000  mm  and  Az  =  0.0625  mm  for  the  air 


Table  3 

Parameters  for  activation  overpotential.  Because  y  in  Aguiar  et  al.  [13]  is  given  as  a 
function  of  temperature,  the  values  calculated  at  7=  750  °C  are  displayed  here 


£act  (kj  mol-1 

) 

y  (xlO9  Am~2' 

) 

Cathode 

Anode 

Cathode 

Anode 

Present  study 

140 

120 

5.0 

5.0 

Campanari  et  al.  [9] 

120 

110 

7.0 

7.0 

Aguiar  et  al.  [13] 

137 

140 

5.2 

28.8 

Costamagna  et  al.  [22] 

160 

140 

7.0 

0.057 
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channel  and  Ax  =  1.5625  mm,  Ay  =  0.1875  mm  and  A z  =  0.0625  mm 
for  the  anode  including  the  fuel  channel.  The  total  number  of  com¬ 
putational  cells  varies  according  to  the  length  of  each  components, 
i.e.,  (y.an  and  /z  rib,  while  the  size  of  each  computational  cell  is  kept 
constant. 

The  local  current  density,  i,  and  the  ohmic  losses,  rjohm,  that  sat¬ 
isfy  Ohm’s  and  Kirchhoff  s  laws  with  the  constitutive  Eqs.  (11  )-(16), 
are  computed  by  using  the  Newton-Raphson  method.  The  electri¬ 
cal  circuit  is  segmented  on  the  same  mesh  as  that  for  the  flow 
simulation  (i.e.,  Ax  =  1.5625  mm  and  Az  =  0.0625  mm).  Inside  the 
electrolyte  and  interconnector,  electric  current  in  the  spanwise 
direction  is  neglected. 

As  mentioned  above,  the  governing  equations  for  the  mass, 
momentum,  and  energy  transports  are  coupled  with  those  for  the 
electric  current  via  the  mass  flux  on  the  TPB  (Eq.  (6))  and  the  heat 
production  (Q.  in  Eq.  (7)  and  q  in  Eq.  (9)).  Iterative  computation  is 
continued  between  the  mass,  momentum  and  energy  transports 
and  the  electric  current  until  the  steady-state  solution  converges 
for  all  the  variables. 

2.4.  Validation 

In  order  to  validate  the  numerical  procedure  described  above, 
we  compare  the  cell  performance  predicted  by  the  simulation  with 
available  experimental  data.  In  most  experiments,  however,  the  sin¬ 
gle  cell  performance  is  evaluated  in  an  isothermal  furnace  and  the 
amount  of  heat  production  is  negligibly  small.  Therefore,  unlike  all 
other  cases  in  the  present  study,  we  neglect  the  heat  production 
term  in  Eqs.  (7)  and  (9),  but  instead  assume  a  uniform  temperature 
of  750  °C.  The  total  length  of  the  cell  is  100  mm  and  both  height  and 
width  of  the  fuel  channel  (ly >ch  and  lz >ch)  is  1.0  mm.  The  thicknesses 
of  the  cathode  and  the  electrolyte  are  100  and  20  |xm,  respectively. 
The  parameters  for  the  cell  structure  are  assumed  as  lz  rib,  =  1.5  mm, 
ly, an  =  1.0  mm,  san  =  0.5  and  dP)an  =  1.0  [xm. 

Fig.  3  shows  the  relationship  between  the  current  and  the  output 
voltage  (/-V)  and  also  between  the  current  and  the  power  density 
(J-P)  obtained  by  the  present  simulation.  A  fairly  good  agreement 
is  found  with  the  experimental  data  [25]  when  the  current  density 
is  relatively  small,  albeit  there  is  some  difference  in  the  configura¬ 
tion  and  operating  condition.  Larger  discrepancy  is  found  for  larger 
current  density,  i.  This  discrepancy  might  have  been  caused  by  the 
strong  sensitivity  of  computational  result  to  the  model  parameters 
that  cannot  be  accurately  determined  from  the  experiments.  For 
instance,  both  the  l-V and  I-P  curves  nearly  collapse  with  the  exper¬ 
imental  data  when  Eact,an  and  the  electrolyte  thickness  /y  eie  are 
adjusted  to  130  kj  mol-1  and  30  [xm,  respectively.  Another  possible 
cause  is  the  applicability  of  the  thin  electrode  assumption,  which 
tends  to  underestimate  anodic  overpotential  when  ionic  resistance 
becomes  dominant  over  activation  overpotential  inside  the  reactive 
porous  electrode  of  finite  thickness  [26]. 

The  grid  spacing  described  above  was  employed  in  this  valida¬ 
tion,  and  it  corresponds  to  the  number  of  computational  cells  of 
64  x  21  x  16  (in  x,  y  and  z  directions)  for  the  air-side  including  the 
cathode,  and  of  64  x  15  x  16  for  the  fuel-side  including  the  anode. 
Simulation  was  also  performed  on  a  computational  grid  twice  finer 
in  each  direction  and  the  resulting  change  in  VceU  was  within  0.5% 
difference. 

3.  Characteristics  of  transport  and  reaction  phenomena 

3.1.  Heat  and  mass  transfer 

The  average  current  density  of  the  electrolyte,  ieleavr,  is 
0.50  A  cm-2.  The  parameters  for  the  anode  structure  are  set  to  be 
/z  rib,  =  1.5  mm,  /y>an  =  1.0  mm,  and  £an  =  0.5  and  dp>an  =  1.0  |xm.  Here- 


Fig.  3.  l-V  and  I-P  characteristics  of  the  isothermal  cell  (7=  750 °C  (constant), 
^gas.in  =  70  atm). 

after,  this  case  is  referred  to  as  the  reference  case.  The  other 
common  parameters  used  throughout  this  study  are  listed  in 
Table  2. 

Fig.  4  shows  the  distributions  of  the  O2  concentration  on  the 
cathode-side  TPB,  H2  and  H20  concentrations  on  the  anode-side 
TPB,  and  electrolyte  temperature,  Teie-  Each  reactant  has  a  concen¬ 
tration  gradient  mainly  in  the  streamwise  (x)  direction,  and  this  is 
characterized  by  the  rates  of  consumption  and  production  of  each 
chemical  species  through  the  reaction.  In  addition,  concentration 
gradients  of  H2  and  H20  are  found  also  in  the  spanwise  (z)  direction, 
especially  on  the  rib  region  of  the  porous  anode. 

Fig.  5  shows  the  normalized  velocity  vector  and  the  molar  frac¬ 
tions  of  H2  and  H20  in  a  cross-stream  (y-z)  plane  at  x  =  50mm. 
The  reactants  are  transported  by  the  cross-stream  flow  induced 
by  the  migration  of  O2-  ions  and  the  diffusion.  The  magnitude  of 
this  cross-stream  flow  is  O  (10-2)  as  compared  to  that  of  the  main 
stream  of  fuel.  The  Peclet  number  of  mass  transfer,  Pem  =  UL/Di  e ff, 
is  O  (10-1),  where  U  and  L  are  the  maximum  velocity  of  the 
cross-stream  flow  and  the  anode  thickness,  ly>an,  respectively.  This 
indicates  that  the  contribution  of  the  cross-stream  flow  to  the  mass 
transfer  in  the  cell  is  negligibly  small  compared  with  that  of  diffu¬ 
sion  in  the  porous  electrodes. 

The  diffusion  paths  of  the  reactants  from  the  air  and  fuel  chan¬ 
nels  to  the  TPB  are  formed  in  the  porous  electrodes.  Eq.  (5)  indicates 
that  the  magnitude  of  reactant  diffusion  inside  the  electrodes  is 
found  to  be  about  ten  times  smaller  than  that  in  the  channel  region. 
As  a  result,  variations  of  H2  and  H20  concentration  are  observed  in 
the  anode  cross-section,  which  lead  the  concentration  gradient  on 
the  TPB.  However,  the  reactant  distribution  on  TPB  in  the  present 
anode-supported  FT-SOFC  is  less  inhomogeneous  than  that  in  the 
conventional  anode-supported  planar  SOFC  [6].  This  is  because  the 
FT-SOFC  has  the  supplementary  diffusion  paths  of  the  reactants 
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Fig.  4.  Molar  fraction  of  02,  H2  and  H20  on  the  TPB  and  electrolyte  temperature  (the  reference  case,  rgasin  =  600  °C,  iele,avr =  0.5  Acrrr2).  The  gas  inlet  is  located  at  x  =  0. 
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(i.e.,  the  rib  of  the  fuel  channel),  which  do  not  exist  in  the  planar 
SOFC. 

Since  the  heat  generated  by  the  exothermic  reaction  is  con¬ 
vened  by  the  air  and  fuel  flows,  Tele  increases  monotonously  from 


the  cell  inlet  toward  the  outlet.  The  temperature  gradient  in  the 
cross-sections  is  small.  The  temperature  distribution  is  determined 
mainly  by  the  heat  conduction  and  the  local  heat  production, 
since  the  Biot  number,  Bi  =  hfUeiJy,anAam  is  considerably  small  (less 
than  0.1  based  on  the  anode  thickness).  The  contribution  of  the 
cross-stream  flow  to  the  heat  transfer  in  the  cell  is  also  small 
as  to  the  mass  transfer,  since  the  Peclet  number  of  heat  transfer, 
PCh  =  PgasCp.gasfff'Agas*  is  O  (10  2). 


3.2.  Electrochemical  reaction 

Fig.  6  shows  the  distributions  of  the  electromotive  force  (EMF), 
Etpb,  and  the  current  density  in  the  electrolyte,  iele.  The  ETPB  takes 
a  maximum  at  the  inlet  and  gradually  decreases  toward  the  cell 
outlet  as  02  and  FI2  are  consumed  and  H20  is  produced.  In  contrast, 
the  location  of  maximum  iele  appears  in  the  middle  of  the  cell.  This 
is  because  the  overpotentials  are  very  temperature  sensitive  and 
there  is  a  trade-off  between  the  increase  of  temperature  and  the 
decrease  of  reactant  concentration.  In  addition,  both  ETPB  and  ie le 
tend  to  be  smaller  in  the  region  above  the  rib.  These  characteristics 
are  similar  to  those  of  the  conventional  planar  SOFCs  with  co-flow 
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Fig.  6.  Electromotive  force  (EMF)  ETPB  and  electrolyte  current  density  iele  (the  reference  case). 
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Fig.  7.  Distribution  of  £bulk  and  spanwise-averaged  overpotentials  (the  reference 
case). 

arrangement,  although  the  spanwise  variations  of  ETPB  and  zeie  are 
found  much  smaller  [10]. 

The  decrease  of  Etpb  and  ie\e  in  the  region  above  the  rib  is  mainly 
due  to  insufficient  diffusion  of  H2  and  H20  in  the  porous  anode. 
However,  as  mentioned  in  Section  3.1 ,  the  H2  shortage  and  the  H20 
stagnation  in  the  anode-supported  FT-SOFC  can  be  alleviated  by 
mass  transfer  through  the  rib  and  the  bottom  of  the  fuel  channel. 
As  a  result,  the  concentration  overpotential  and  the  spanwise  varia¬ 
tions  of  Etpb  and  ie le  become  smaller  than  those  of  the  conventional 
planar  SOFCs.  These  characteristics  lead  to  an  advantage  of  FT-SOFC 
over  planar  SOFC  under  the  conditions  of  high  fuel  utilization. 

The  componential  losses  along  the  streamwise  direction  are 
plotted  in  Fig.  7.  It  reveals  that  the  dominant  components  for  the 
voltage  loss  are  the  activation  overpotential  (rj act)  at  the  electrodes, 


the  concentration  overpotential  (rjc on),  and  the  ohmic  loss  of  the 
electrolyte  (^ohm.eie)-  The  ohmic  losses  of  the  electrodes  {rjohmca 
and  ?7ohm,an)  and  the  interconnector  (^0hm,int)  are  negligibly  small 
due  to  their  small  resistivities.  The  cathodic  activation  overpotential 
(^act.ca)  and  ?7ohm,eie  are  maximized  at  the  cell  inlet.  This  is  natural 
because  these  losses  strongly  depend  on  the  temperature  and  sim¬ 
ply  decrease  as  the  temperature  increases.  In  contrast,  the  anodic 
activation  overpotential  (?7act,an),  which  should  be  determined  by 
the  combination  of  zeie,  the  local  concentrations  of  H2  and  H20,  and 
the  temperature,  remains  almost  constant  over  the  streamwise  cell 
length. 

3.3.  Parametric  survey 

3.3.1  Effect  of  the  average  current  density 

We  examine  the  condition  at  which  the  average  current  density 
is  a  half  of  the  reference  case,  i.e.,  zeie,avr*  =  0.25  A  cm-2.  The  other 
parameters  remain  unchanged  from  the  reference  case.  The  result¬ 
ing  output  voltage,  Vcen ,  is  0.769  V  (cf.,  Vcen  =  0.661  V  in  the  reference 
case). 

Fig.  8  shows  the  distributions  of  the  electrolyte  temperature, 
Teie,  and  the  local  current  density,  iele.  It  is  found  that  the  lower 
z’eie.avr  makes  the  temperature  higher  near  the  inlet  and  lower  near 
the  outlet.  The  location  of  maxinum  zele  shifts  toward  upstream 
with  the  decrease  of  ieie,avr-  The  lower  temperature  near  the  outlet 
is  simply  due  to  the  smaller  amount  of  total  heat  generation.  The 
higher  temperature  near  the  inlet  is  caused  by  deterioration  of  heat 
exchange,  because  the  air/fuel  flow  rate  changes  proportionally  to 
ieie.avr  under  the  condition  of  constant  air/fuel  utilization  rate. 

The  concentrations  of  H2  and  H20  in  the  y-z  plane  at  x  =  50  mm 
are  shown  in  Fig.  9.  Due  to  the  reduced  local  reaction  rate,  the 
concentration  gradient  in  the  cross-stream  plane  is  mitigated  as 
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Fig.  8.  Electrolyte  temperature  and  current  density  of  the  electrolyte  with  ieie,avr =  0.25  A  cm-2. 
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Fig.  10.  Distribution  of  Ebuik  and  spanwise-averaged  overpotentials  at 

iele,avr  =  0-25Acm-2. 

compared  to  the  reference  case.  This  enhancement  of  fuel  supply 
leads  to  the  smaller  spanwise  variation  of  iele,  as  is  observed  in 
Fig.  8,  and  the  smaller  concentration  gradient  decreases  both  r]con 
and  ?7act,an  as  shown  in  Fig.  10.  Moreover,  the  lower  ieie,avr  and  the 
higher  temperature  near  the  inlet  result  in  the  decrease  of  pac t,  both 
for  the  anode  and  the  cathode,  and  also  of  ^ohm.eie-  As  a  combina¬ 
tion  of  these  effects,  the  output  voltage  increases  with  the  decrease 
of  ieie 

,avr- 

3.3.2.  Effect  of  the  design  parameter 

Several  cases  with  different  design  parameters  of  the  anode,  i.e., 
the  pore  diameter,  dp>an,  the  porosity,  £an,  the  rib  width,  Zzrib,  and 
the  rib  thickness,  (y>an,  are  studied.  Table  4  summarizes  the  resulting 
output  voltage,  Vcen,  for  these  cases.  As  compared  to  the  reference 
case,  VceU  is  significantly  improved  in  the  case  of  larger  dp,an,  while 
only  marginally  higher  in  the  other  cases  studied. 

The  increase  of  Vcell  is  explained  by  the  mass  transfer  enhance¬ 
ment  in  the  cross-sectional  plane  as  shown  in  Fig.  11.  In  all  cases, 
the  fuel  concentration  in  the  porous  anode  becomes  uniform  and 
that  on  the  TPB  is  larger  than  the  reference  case.  This  enhanced 
diffusion  improves  the  supply  of  FI2  and  mitigates  the  stagna¬ 
tion  of  H20,  so  that  rjcon  and  ?7act,an  are  reduced.  In  the  cases  of 
larger  dp,an  or  £an,  the  enhancement  of  diffusion  is  attributed  to  the 
increase  of  the  effective  diffusivity,  Di  e ff.  As  suggested  by  Eq.  (5), 
Di  eff  increases  with  increase  of  dp>an,  because  the  Knudsen  diffu¬ 
sivity,  Dj  k,  is  proportional  to  dp>an.  Eq.  (5)  also  indicates  that  Dieff 
increases  as  £an  is  increased.  On  the  other  hand,  for  the  anode 
with  a  thinner  or  narrower  rib,  shorter  diffusion  length  directly 
results  in  the  enhancement  of  mass  transfer,  although  Di  eff  is 
unchanged. 

Fig.  12  shows  the  distribution  of  the  electrolyte  temperature.  The 
boundary  and  operating  conditions  are  to  the  same  as  those  of  refer¬ 


ence  case.  The  temperature  drops  with  the  increase  of  £an,  whereas 
it  is  almost  unchanged  by  the  variation  of  dp>an.  As  discussed  in 
Section  3.1,  the  heat  transfer  in  the  porous  anode  is  dominated  by 
the  streamwise  heat  conduction.  It  is  nearly  proportional  to  the 
cross-sectional  area,  San,  and  the  effective  thermal  conductivity, 
Xeff,  which  is  defined  as 

A-eff  =  £Agas  +  (1  —  s)XS0\.  (17) 

Therefore,  the  observation  in  each  case  can  be  explained  as  fol¬ 
lows: 

•  For  larger  pore  (dp,an  =  2.0  fxm): 

The  temperature  distribution  is  unchanged  because  San  and  Aeff 
are  independent  of  dp,an. 

•  For  higher  porosity  (£an  =  0.6): 

The  heat  conduction  is  reduced  because  Agas  is  much  smaller 
than  Xso\.  The  reduced  heat  conduction  causes  the  temperature 
drop  and  the  increase  of  pac t  and  rj 0hm,eie  near  the  inlet.  The  deterio¬ 
ration  due  to  this  temperature  drop  is  nearly  canceled  by  the  favor¬ 
able  effect  due  to  the  mass  transfer  enhancement  discussed  above. 

Fig.  12  also  shows  that  the  inlet  temperature  decreases  with  the 
decrease  of  lz  rib  and  Zy,an.  When  the  rib  is  thin  or  narrow,  the  stream- 
wise  heat  conduction  is  suppressed  due  to  the  reduction  of  the 
cross-sectional  area.  This  leads  to  the  decrease  of  inlet  temperature 
and  the  increase  of  pact  and  ^0hm,eie»  similarly  to  the  case  of  large 
£an.  In  these  cases,  the  decrement  of  cross-sectional  area  is  compa¬ 
rable  to  that  of  diffusion  length.  Thus,  a  trade-off  occurs  between 
the  mass  transfer  enhancement  and  the  heat  conduction  degrada¬ 
tion  when  Zzrib  and  Zy,an  are  varied.  As  a  result,  Vcell  is  relatively 
insensitive  to  these  parameters. 

In  the  present  study,  the  activation  overpotential  model,  Eq.  (15) 
is  fixed  during  the  parametric  survey  of  the  anode  microstructure, 
i.e.,  dp,an,  and  £an.  Flowever,  recent  efforts  have  shown  that  intro¬ 
ducing  thin  anode  functional  layer  can  improve  both  activation  and 
concentration  overpotentials  [27-29].  Besides,  it  is  estimated  that 
the  effective  anode  reactive  thickness  for  standard  grain  size  (1  p,m) 
is  around  10  p,m  [30],  which  is  much  smaller  than  the  total  anode 
thickness.  Thus,  it  seems  reasonable  to  assume  thin  active  layer 
and  to  optimize  bulk  anode  support  layer  independently  without 
changing  the  activation  overpotential  at  the  anode  electrolyte  inter¬ 
face. 

4.  Development  of  simplified  cell  model 

4.1.  Model  description 

In  this  section,  a  simplified  numerical  model  for  the  anode- 
supported  FT-SOFC  is  developed  based  on  the  knowledge  obtained 
by  the  3D  simulation  presented  above.  The  results  of  3D  simula- 


Table  4 

Dependency  of  Vcen  on  ieje  avr  and  anode  design  parameters  obtained  by  3D  and  simplified  simulations 


Values  in  bold  are  differences  from  those  of  the  reference  case. 
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Fig.  11.  Normalized  velocity  vectors  and  molar  fraction  of  H2  in  a  cross-section  x  = 
anode;  (d)  narrower  rib.  See,  Table  4  for  the  specification  of  each  case. 


Fig.  12.  Electrolyte  temperature  under  different  design  parameters 
Oele.avr  =  0.5  A  cm  ^ ). 


tion  confirm  that  lateral  (z)  variation  of  the  reaction  field  at  TPB 
for  the  anode-supported  FT-SOFC  is  very  small,  and  the  advec- 
tion  terms  of  the  transport  equations  can  be  neglected.  Namely, 
the  distribution  of  the  electrochemical  field  (i.e.,  current  density, 
EMF,  and  overpotentials)  at  TPB  can  be  assumed  to  be  uniform  in 
the  spanwise  direction,  and  the  air  and  fuel  streams  to  be  one¬ 
dimensional  plug  flows.  The  resulting  model  is  basically  similar  to 
the  conventional  one-dimensional  model  [11,13,15,31,32].  In  this 
model,  a  special  procedure  is  adopted  to  treat  the  concentrations 
of  H2  and  H20  in  the  thick  anode  with  a  complex  cross-sectional 
geometry. 

The  cell  is  divided  into  six  segments,  as  shown  in  Fig.  13(a).  These 
segments  as  well  as  the  air  and  fuel  channels  are  discretized  by 
control  volumes  in  the  streamwise  (x)  direction.  The  physical  quan¬ 
tities,  i.e.,  the  mass  flux  of  each  species  and  the  temperature,  are 
defined  at  each  control  volume.  The  solution  for  the  set  of  conser¬ 
vation  equations  (described  in  the  following  sections)  is  obtained 
by  using  the  Newton-Raphson  method. 
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Fig.  13.  (a)  Segmented  component  model  of  the  anode-supported  FT-SOFC.  (b)  Distribution  of  0  in  the  cross-sectional  area  of  the  anode  and  fuel  channel  (Reference  case). 
The  thin  layers  of  electrolyte,  cathode,  and  interconnector  (of  which  thicknesses  are  20  p,m,  100  p,m  and  100  p,m,  respectively,  and  treated  as  boundaries  in  the  simulation) 
are  drawn  not  in  scale  so  as  to  be  visible. 


4.1.1.  Transport  in  fuel  and  air  channels 

The  transports  of  the  mass-flux  of  each  reactant,  rhj,  and  the 
temperature,  T,  in  the  fuel  and  air  channels  are  expressed  as 


J  -^=±('z,ch  +  (z,rib)jj  j  S02,  H2,H20 

I  9/71/ 

(-g^=0  jeN2,  C02 

and 


(18) 


_9 

dx 


dx  —  C/CC)nd  +  Qc  onv,x  +  Qconv,y  +  Qht 


air  side:  jeN2,02  ke  02 

fuel  side:  jeH2,H20,C02  ke  H2,H20, 


(19) 


where 


4cond  —  A:h 


dx, 


(20) 


dx, 

T 

Cp,lJ, 


Qht  —  ^ht^C^gas  Tsoi). 


(21) 


(22) 

(23) 


Here,  Ach  and  Aht  are  the  cross-sectional  area  of  the  flow  channel 
and  the  area  that  participate  in  the  convective  heat  transfer,  respec¬ 
tively.  Each  term  on  the  right  hand  side  of  Eq.  ( 19)  expresses  the  heat 
conduction,  the  convective  heat  transfer,  and  the  heat  convection 
on  the  TPB.  In  order  to  determine  the  heat  transfer  coefficient,  h, 
in  Eq.  (23),  the  Nusselt  number  for  the  parallel  plates  (NuH  =  8.235) 
and  that  for  the  square  duct  (NuH  =  3.599)  [33]  are  used  for  the  air 
and  fuel  channels,  respectively. 


4.1.2.  Transport  in  solid  and  porous  elements 

As  mentioned  in  the  Section  3.1 ,  the  Peclet  number  of  heat  trans¬ 
fer  in  each  cross-sectional  area  of  the  cell  is  considerably  small. 
In  addition,  the  main  streams  of  air  and  fuel  hardly  affect  the 
heat  transfer  inside  the  electrodes.  Namely,  the  heat  transfer  in 
the  electrodes  is  governed  by  the  heat  conduction  only.  With  this 


assumption,  the  heat  transport  in  the  element  j  (see,  Fig.  13(a))  can 
be  expressed  as 


\YCJ 


9  ( ■.  9TS0|j 
dx  {Xef(J~dT 


dx  —  Qhtj  +  )  Qcond.j.k  +  Qgene.j ' 


(24) 


where 

„  _  f  4htjh(7gas  —  TS0]j)  j  e  1 , 3,  4,  5,  6 

%tJ  “jo  je  2 


(25) 


ijAy,j  ^ohm,j 


AH, 


H20 


2  F 


I  +  (^ohmj  +  ^act) 


j  g  1, 3,  4,  5,  6 
je  2 


ny,j,k 


dL 


sol 


(Sj/Xeffj)  +  (9k/Aeff?k)  9y 


(26) 


(27) 


Here,  Axj  and  Ayj  are  the  areas  of  heat  transmission  in  the  stream- 
wise  and  perpendicular  directions,  respectively,  and  8j  is  the 
thickness  of  the  component  j.  In  Eq.  (24),  the  electrodes  are  treated 
as  solid  with  an  effective  heat  conductivity  of  A.eff.  The  summation 
of  Eq.  (25)  for  j  is  equivalent  to  Eq.  (23).  Note  that  the  effect  of  the 
anode  design  can  be  investigated  by  changing  the  parameters  of 
microstructure,  £an  and  dp,an,  and  the  size  of  elements  3,  4  and  5 
shown  in  Fig.  13(a). 


4.1.3.  Concentration  of  reactant  on  TPB 

In  the  conventional  models,  the  molar  concentration  that  par¬ 
ticipates  in  the  electrochemical  process  at  the  TPB  is  given  by 


Xi, tpb  =  ^j,ch  ± 


*ele  KqT  1 

nF  PD i  eff  Zd jff 


(28) 


where  n  is  the  number  of  electrons  in  the  reaction,  while  Zdiff  is  the 
diffusion  length  being  equal  to  the  anode  thickness  in  the  case  of 
planar  SOFC  [34].  However,  Zdiff  cannot  be  determined  for  arbitrary 
cross-sectional  geometry.  Also,  as  observed  in  Figs.  5, 9  and  11 ,  vari¬ 
ations  of  H2  and  H20  concentrations  in  the  cross-section  of  anode 
cannot  be  neglected.  Hence,  the  following  procedure  is  employed 
to  obtain  the  spanwise-averaged  concentration  of  reactants  on  the 
TPB. 

We  assume  that  the  physical  properties  are  constant  in  the 
cross-section  and  that  the  concentration  variation  in  the  stream- 
wise  direction  is  much  milder  than  that  in  the  cross-section.  Then, 
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the  transport  equation  for  the  species,  Eq.  (4),  reduces  to  the  two- 
dimensional  Laplace  equation  by  neglecting  the  advection  term, 
i.e., 


9 y2  9z2 

Here,  the  normalized  mass  fraction,  cp,  is  defined  as 


(29) 


0(y,  z) 


yj(y,z)-y,jPB 
yi,ch  -  yj.TPB 


(30) 


with  Yl)Ch  and  YIiTPB  being  the  mass  fractions  of  reactant  i  in  the 
flow  channel  and  that  on  the  TPB,  respectively.  This  Laplace  equa¬ 
tion  is  solved  with  boundary  conditions  of  0  =  1  in  the  fuel  channel 
and  0  =  0  on  the  TPB.  The  spanwise  (z-)  gradient  of  0  is  zero  on 
the  symmetric  planes  (i.e.,  for  anode-side,  90/9z  =  0  at  z  =  0  mm  and 
z  =  2  mm  in  Fig.  13(b)).  Once  Eq.  (29)  is  solved  for  0,  as  exemplified 
in  Fig.  13(b),  the  gradient  of  Yz  at  the  TPB  can  be  computed  as 


dYi 

dy 


=  (^i,ch  -  ^i,TPb) 

TPB 


90 

dy 


(31) 


With  an  assumption  of  constant  physical  property,  the  molar 
fraction  of  reactant  i  on  the  TPB  can  be  computed  by  using  Eqs.  (6), 
(12),  (31)  and  the  gas  state  equation,  as 


Xi,  tpb  =  ^i,ch  =•= 


*ele  RoT 

nF  PD i  eff 


(32) 


The  normalized  concentration  gradient  on  the  TPB  (90/9y) |Tpb* 
depends  only  on  the  cross-sectional  geometry  of  the  cell  and  it 
can  be  obtained  by  solving  the  Laplace  equation  only  once  for  the 
given  geometry.  Note  that  Eq.  (32)  is  equivalent  to  Eq.  (28)  for  the 
cathode-side,  if  l/(90/9y)TPB  is  replaced  by  lytC a. 


4 A. 4.  Electrochemical  reaction 

The  governing  equations  and  numerical  scheme  for  the  electro¬ 
chemical  reaction  is  the  same  as  in  the  case  of  3D  simulation.  The 
only  difference  in  the  simplified  model  is  the  spanwise-averaged 
equivalent  circuit  model,  which  has  a  single  current  path  in  a  set  of 
electrolyte  and  electrodes. 


case). 


4.2.  Validation  of  the  model 

Fig.  14  shows  the  distributions  of  ie ie  and  Tele  predicted  by  the 
simplified  model  in  the  cases  of  ieie,avr  =  0-25  and  0.5  A  cm-2  (refer¬ 
ence  case).  The  profiles  computed  by  using  the  simplified  model  are 


Rib  width  lzrib  [mml  Rib  width  lz  rib  [mm]  Rib  width  lz  rih  [mm] 


Fig.  15.  Output  voltage  (upper  row)  and  volumetric  power  density  (lower  row)  for  wide  range  of  parameters:  (a)  ieieavr  =  0.25  Acm  2;  (b)  iejeavr  =  0.50Acm  2;  (c) 
lele,avr  =  0.75  A  cm-2.  The  open  circle  in  (b)  corresponds  to  the  reference  case. 
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in  good  agreement  with  those  of  3D  simulations.  In  the  reference 
case,  the  deviations  from  the  3D  simulation  are  within  9  mA  cm-2 
for  the  maximum  current  density,  10.4  K  for  the  maximum  tem¬ 
perature,  and  within  5  mm  for  the  peak  location  of  current  density. 
This  difference  may  be  due  to  the  assumptions  of  constant  Nusselt 
number  and  spanwise  uniformity  of  the  reaction  rate.  The  Nus¬ 
selt  number  is  much  larger  in  the  thermal-entry  region  [35],  of 
which  length  is  increased  with  increase  of  ieie,avr-  The  assumption  of 
constant  Nusselt  number  results  in  overprediction  of  the  tempera¬ 
ture  near  the  inlet  and  underprediction  of  the  overpotential.  On  the 
other  hand,  as  discussed  in  Section  3.3,  the  spanwise  variation  of  zele 
is  pronounced  when  the  mass  diffusion  is  insufficient.  Namely,  the 
assumption  of  uniform  reaction  rate  might  have  caused  the  errors 
in  ?7con»  Pact  ^nd  ieje. 

The  resulting  output  voltage,  \7cen,  is  shown  in  Table  4  for  dif¬ 
ferent  values  parameters.  Although  the  computed  voltage  is  in 
fair  agreement  with  the  3D  simulation,  discrepancy  is  also  found 
depending  on  the  parameters.  Agreement  with  the  3D  model  is 
found  to  be  better  for  cells  with  a  thick  anode  with  narrow  rib,  in 
which  the  spanwise  variation  of  reaction  is  relatively  small.  The 
computational  cost  of  this  simplified  model  is  1/105  that  of  the  3D 
simulation.  These  results  suggest  that  present  simplified  model  of 
the  anode-supported  FT-SOFC  can  be  used  as  a  convenient  tool  for 
the  parametric  survey  of  cell  design  and  operating  conditions. 

4.3.  Example  application 

As  an  example  application  of  the  simplified  model,  we  perform  a 
parametric  survey  with  wide  ranges  of  anode  thickness,  ly<an  (from 
0.2  to  3.0  mm  with  increment  of  0.2  mm),  rib  width,  Zzrib  (0.25  to 
4.00  mm  with  increment  of  0.25  mm),  and  average  current  density, 
ze le.avr  (0.25, 0.5  and  0.75  A  cm-2 ).  All  other  parameters  are  assumed 
to  be  the  same  as  those  of  the  reference  case. 

Fig.  15  shows  the  computed  voltage,  Vcen,  and  the  volumetric 
power  density,  Wcell,  defined  as 

Wce„  =  Vce,"lele'avr,  (33) 

Zy,cell 

where  Zycell  is  the  cell  thickness,  which  is  the  summation  of 
electrolyte,  electrodes  and  interconnector  thicknesses.  As  can  be 
noticed  from  this  figure,  we  can  find  a  set  of  parameters  to 
attain  desired  performance  by  using  the  present  simplified  model: 
for  instance  (ly> an,  lz, rib) =  0-4,  0.5  mm)  for  the  maximum  Vcell  at 
z’eie.avr =  0-5  A  cm-2,  and  (ly> an,  ^z.rib) =  (0.2, 0.8  mm)  for  the  maximum 
Wceu  ^  Zeie.avr  =  0.5  A  Cm~2. 

5.  Conclusions 

Numerical  simulation  and  modeling  of  the  heat/mass  transfer 
and  electrochemical  reaction  in  an  anode-supported  flat-tube  solid 
oxide  fuel  cell  (FT-SOFC)  are  carried  out.  Through  the  numerical 
analyses  with  a  parametric  survey,  the  following  conclusions  are 
derived: 

1.  The  distributions  of  the  reaction  fields  in  the  anode-supported 
FT-SOFC  are  quantitatively  similar  to  those  of  the  conventional 
planar  cell  with  co-flow  arrangement.  However,  the  concentra¬ 
tion  and  activation  overpotentials  are  mitigated  in  FT-SOFC  by 
enhanced  reactant  diffusion  in  the  porous  rib  of  the  fuel  channel. 
These  characteristics  give  an  advantage  to  FT-SOFC  over  planar 
SOFC  under  the  condition  of  high  fuel  utilization  and  high  current 
density. 

2.  The  design  parameters  affect  the  performance  characteristics 
as  well  as  flow  field  of  the  cell.  For  a  fixed  activation  over¬ 
potential  model,  a  higher  output  voltage  is  obtained  by  using 


an  anode  with  larger  pore  diameter.  However,  the  output  volt¬ 
age  of  the  cell  is  less  sensitive  to  the  cross-stream  geometry  of 
anode  due  to  a  trade-off  between  the  heat  transfer  and  the  mass 
transfer. 

3.  A  simplified  numerical  model  of  the  anode-supported  FT-SOFC 
is  developed  based  on  the  results  obtained  by  the  three- 
dimensional  simulation.  This  model  can  predict  sufficiently  well 
the  output  voltage  as  well  as  the  distributions  of  temperature 
and  current  density  with  very  low  computational  cost.  Thus,  it 
can  be  used  as  a  convenient  tool  for  surveying  wide  range  of 
anode  support  FT-SOFC  design  parameters.  In  order  to  determine 
the  optimal  design  of  the  anode-supported  FT-SOFC,  however, 
further  parametric  survey  should  be  required. 
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